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Thermalisation of configurations with initial white noise power spectrum is studied in numerical 
simulations of a classical one-component $ 4 theory in 2+1 dimensions, coupled to a small amplitude 
homogenous external field. The study is performed for energy densities corresponding to the broken 
symmetry phase of the system in equilibrium. The effective equation of the order parameter motion 
, is reconstructed from its trajectory which starts from an initial value near the metastable point and 

■ ends in the stable ground state. This phenomenological theory quantitatively accounts for the decay 

of the false vacuum. The large amplitude transition of the order parameter between the two minima 
displays characteristics reflecting dynamical aspects of the Maxwell construction. 
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I. INTRODUCTION 

> ■ 

. The reach of equilibrium from a metastable state covers a large number of interesting effects from instabilities 
' observed in the mixed phase of first order phase transitions JIJ to the inflation in the early Universe The final 
, state is reached in an irreversible process, the description of this relaxation relies on out-of-equilibrium dynamics of 
■^sl" ' field theories. 

The decay of metastable states is usually discussed in the framework of the nucleation scenario ||. It has been 
implemented in the form of saddle point expansions in classical [Q, and quantum systems J^]. This large amplitude 
instability is the first of the possible instabilities, suggested by mean field analysis, which consists of the creation of 
a bubble of the true vacuum larger than the critical size, embedded into the false one. 

Another possibility, the instability against fluctuations with infinitesimal amplitude leads to the spinodal phase 
O-i separation. A recent observation made it clear that soft fluctuations of these inhomogeneous unstable modes gen- 
£h ' erate in equilibrium the Maxwell construction by their tree-level contribution to the renormalization group flow jfj. 
The fluctuation induced flatness of the effective potential suggests the dominance of spinodal phase separation in 
equilibrium. 

Since the type of instability one observes, might depend essentially on the time scale of the observation, a detailed 
investigation of the time evolution can separate the effects of the two kinds of instabilities. This is made possible by 
C3 large scale computer simulations of the thermalisation process in closed systems. 

Another question, left open by the Maxwell construction in equilibrium |?]], concerns the structure of the vacuum 
with spontaneous symmetry breaking. In fact, the effective potential is related to the probability distribution of 
the order parameter and the Maxwell-cut applied to the former suggests that the latter is also degenerate in the 
mixed phase. Either we accept that the vacua with spontaneously broken symmetry correspond overwhelmingly to 
the mixed phase or a dynamical mechanism is seeked to eliminate the mixed phase from among the final states of the 
time evolution. 

In cosmology different slow-roll scenarios of inflation are being considered. Recent studies of the dynamics of inflaton 
fields with large number of components (large N limit) displayed for a first time a dynamical version of the Maxwell 
construction J8| . The Hartree type solution of the quantum dynamical equations leads to the conclusion that the order 
parameter might get rest with finite probability at any value smaller than the position of the stable minimum of the 
tree level effective potential, corresponding to the stabilisation of a mixed state. 

Detailed investigations of the thermalisation phenomenon were performed also in noisy-damped systems, coupled 
to an external heat bath P-[T^]. The relaxation to thermal equilibrium of the space averaged scalar field (the order 
parameter) starting from metastable initial values has been thoroughly investigated. In these simulations the damping 
coefficient is treated as an external control parameter. Using the numerical solution of the corresponding Langevin- 
equations the validity range of the analytical results for the homogenous nucleation mechanism has been explored. 
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In this paper we focus on an alternative description of the decay process of the metastable vacuum state. The process 
is described exclusively in terms of the homogenous order parameter (OP) mode. The evolution of the OP is studied in 
interaction with the rest of the system as described by the reversible dynamical equations of motion of the full system. 
Careful analysis of its dynamics allows us to explore the effects of both kinds of the above mentioned instabilities. 
The transition of the order parameter from the metastable to the stable vacuum is induced by a homogenous external 
"magnetic" field, whose strength is systematically reduced. No random noise is introduced to represent any external 
heat bath, the friction coefficient of the effective order parameter dynamics is determined internally. 

Our results offer a " dualistic" resolution of the competition between the nucleation and the spinodal phase separation 
mechanisms in establishing the true equilibrium. On the one hand, we find that the statistical features of the decay of 
the false vacuum agree with the results obtained by expanding around the critical bubble. The microscopic mapping 
of the field configurations during the relaxation supports the bubble creation scenario. Alternatively, the effective 
OP-thcory displays the presence of soft modes and produces dynamically a Maxwell-cut when the time dependence of 
the transition trajectory is described in the effective OP theory. We find that the larger is the system the smaller is 
the external field which is able to produce the instability. For infinite systems an infinitesimal field pushes the system 
through the Maxwell-cut, where no force is experienced by the OP. Therefore it will not stop before reaching the true 
homogenous ground state passing by the mixed states with constant velocity. 

Our model, a spacelike lattice regulated classical scalar $ 4 field theory in its broken symmetry phase is introduced in 
Section 2. In Section 3 we describe the evolution of the system starting from order parameter values near a metastable 
point which relaxes first to this state under the combined effect of parametric resonances and spinodal instabilities. 
The second stage of the transition to the stable ground state is the actual focus of our paper. 

Characteristic intervals of the observed order parameter trajectory are reinterpreted as being the solutions of some 
effective point-particle equation of motion, which displays dissipation effects explicitly. Our approach can be under- 
stood also as the real-time version of the lowest mode approximation used for the estimation of finite size dependences 
in Euclidean field theory [|l 3|J 1 -4| 1 . In this sense our approach can be considered also as the numerical implementation 
of a real time renormalisation group strategy. 

One of our principal goals is to reconstruct the thermodynamics of the classical " OP-ensemble" on the (meta)stable 
branches of the OP-trajectory (Section 4). Its dissipative dynamical equations near equilibria will be established. On 
the transition trajectory we shall elaborate on the presence of the Maxwell construction in the effective dynamics 
describing the motion after nucleation (Section 5). The statistical aspects of the approach to the equilibrium are 
established for reference and comparison in an Appendix. The conclusions of this investigation are summarised in 
Section 6. 

The results of this study can be usefully compared with classical investigations of metastability and nucleation in 
the kinetic Ising model |l5| . This system has first order dissipative dynamics by its definition. Still several relaxation 
features of the kinetic Ising model are comparable to our findings, since the "numerical experiment" performed in the 
two models are essentially the same. Especially, the relaxation function of is simply related to the order parameter 
we focus our attention. In both cases in the mechanism of the bubble growth the aggregation of spontanously generated 
small size regions of the true ground state to its surface plays important role. 



The energy functional of a classical system in a two-dimensional box of size Ld coupled to an external magnetic 
field of strength hd is of the following form: 



The index d is introduced to distinguish the dimensionfull quantities from the dimensionless ones, defined by the 
relations (for m? < 0): 



II. CLASSICAL CUT-OFF $ 4 -THEORY ON LATTICE 




(2.1) 



t = t d \m\, 



x = Xd\m\, 




(2.2) 



For the spatial discretisation one introduces a lattice of size 

L d = Na d = Na—. 

\m\ 

The energy functional of the lattice system can be written as 
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(Here we have introduced the dimensionless time-step a t , which should be chosen much smaller than a, and n denotes 
the lattice site vectors.) The equation of motion to be solved numerically is the following: 

+ a t ) + * n (f - a t ) - 2$ n (t) - i £«(*„+!(*) + * n _i(t) - 2*„(t)) 

+ ^(-$ n + $3-/i) = 0. (2.5) 



The initial conditions for Eq.(2.5) were chosen as 

K(t = 0) = 0, $ n (f = 0) = $ + (2.6) 

The random variable £ is distributed evenly on the interval (—1/2, 1/2). Therefore the starting OP-value is <f>o. The 
energy density E/Na 2 is controlled through the magnitude of <£>i. In this study we have chosen <f>o = 0.815 and 
$i = 4/v / 6- The latter corresponds to a temperature value Tj = 0.57 in the metastable equilibrium (from Eq.(jO])). 
This value is much below the critical temperature of the system (T c ~ 1.5Ti). It has been checked that at this energy 
density all other choices of <E>o > 0, for fixed h, find a unique metastable equilibrium. 



Eq.(2__5) was solved with a = 1 and with typical at values in the range (0.01-0.09). It has been checked that 
the statistical characteristics of the time evolution is not sensitive to the variation of a t , though the "release" time 
(the moment of the transition from metastability to the true ground state) in any single run with given initial 
conditions might change considerably under the variation of a t /a. Three lattice sizes were systematically explored: 
N = 64, 128, 256. Several single runs were performed also for N = 512 and N = 1000 with the aim to analyze in 
more detail some self-averaging physical quantities on different portions of the trajectory under the assumption of the 
ergodicity of the system. The magnetic field h inducing the transition was varied in the range h E -(0,0.08)/\/6. For 
the reconstruction of the effective potential felt by the OP also positive values were chosen up to h = 0.5/ VE. The 
smaller the value of \h\ was fixed, the longer the "release" times have grown on the average. For this reason also the 
runs were prolongated with decreasing \h\, and for the smallest \h\ the length of a run reached up to (10 6 — 10 7 )|m| _1 
until the transition took place. 

For a careful comparison of our transition statistics with the generally used statistical nucleation theory, and also 
for understanding the systematics of its change when h has been diminished a large number of (h, N) pairs were 
used in this analysis. Altogether 24 422 transition events have been recorded (for N = 64: 16 908, N = 128: 2 903, 
N = 256: 4 411). For the largest systems at the smallest h the event rate was 1-2/day/ 400 MHz-processor. 



III. TIME-HISTORY OF THE ORDER PARAMETER 



A typical OP-history is displayed in Fig.gj In the same figure we show also the history of the OP mean square 
(MS)-fluctuation {(Q 2 ) — ( C E > ) 2 ) and of its third moment (((<& — ($)) 3 )). The evolution of the non-zero k modes is 
demonstrated in Fig.pl where the averaged kinetic energy content of the |k| < 2.5 and |k| > 2.5 regions is followed. 
Although the separation value is somewhat arbitrary, namely it divides into two nearly equal groups the spatial 
frequencies available in the lattice system, the figure demonstrates the most important features of the evolution of the 
power in the low-|k| and high-|k| modes. 

In general, five qualitatively distinct parts of the trajectory can be distinguished, although some of the first three 
might be missing for some initial configurations and/or magnetic field strengths. 
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The time evolution of the order parameter, its MS fluctuation and the third moment. The example is selected from 
TV = 512 lattice with h — — 0.04/v / 6 external source strength. 



The OP-motion usually starts with large amplitude damped oscillations. The "white noise" initial condition of 
Eq. ( |2.6| ) corresponds to a k- independent Fourier amplitude distribution, therefore the initial distribution of the kinetic 
energy is ~ u> 2 (|k|). During this period, in the power spectrum of the kinetic energy, first a single sharp peak shows up 
at a resonating |k|-valuc (|k| ~ 1.5), which breaks up into several peaks (|k| < 1) at later times due to the non-linear 
interaction of the modes. At the end of the first period the whole |k| < 1 range gets increased power, the |k| > 1.5 
part of the power spectrum does not seem to change. 
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FIG. 2. The time evolution of the kinetic energy content of the |k| > 2.5 and |k| < 2.5 regions averaged over the corresponding 
k|-intervals. The example is the same as in Fig. 1. 

Next a slow, almost linear (modulated) decrease of the OP follows. At the same time its MS fluctuation increases 
linearly. In Fig.|^ an energy flow towards the low-|k| part can be observed, which proceeds through the excitation 
of single modes in this part. As a result the minimum of the effective potential is continously shifted to smaller 

values, as if the temperature would gradually increase (see the left part of Fig.Q). The full (microscopic) kinetic 
energy density shows in this period less than 5% variation. In view of the picture based on the Maxwell construction 
it is rather surprising that independently of the initial conditions the OP converges towards a well-defined absolute 
value, depending only on the total energy density. 

On the third portion the average value of the OP and its moments stay constant (see the right picture in Fig.[j]). The 
average energy content of the low and high-|k| part of the spectra is nearly the same. This suggests the establishment 
of a sort of thermal (meta)cquilibrium. 

In terms of the terminology introduced for the inflation, the first period leading to this (quasi) -stationary state 
can be called preheating, and the second reheating. Directly before the moment of the transition to the true vacuum 
a peak appears in the power spectrum in the narrow neighbourhood of k = with varying position in time. On the 
snapshots of the real space configurations a set of randomly distributed small bubbles of the true ground state appear 
with a radius increasing in time, until one of the bubbles exceeds the critical size. 

The fourth portion of the motion is the transition itself. The value of the OP MS-fluctuation increases by about a 
factor of three and the temporal width of this transient increase measures very well the transition time. The transition 
time decreases on larger lattices, the height of the jump in the OP-fluctuation is not sensitive to the lattice size. Also 
the third reduced moment shows a characteristic variation. An increase of the temperature proceeds smoothly during 
the transition of the OP to its stable value (Fig.||) . The slight separation of the two curves in Figj^ gives a feeling on 
the degree of uniformity of the temperature variation of the different modes. 

The last portion of the trajectory represents stable (thermaliscd) oscillations around the true ground state. Here 
a complete equilibration of the power spectrum can be observed (see FigJ^) corresponding to a somewhat increased 
temperature (T ~ 0.6). 



IV. MOTION NEAR THE (META)STABLE POINT 

Our analysis of the motion around the (meta)stable value of the order parameter explores the consequences of 
assuming the crgodicity hypothesis for a sufficiently long finite time interval, after the system has already reached 
the equilibrium. The equilibrium is characterised by a limiting probability density in the configuration space. The 
averaging with this density should provide the same value as the one yielded by a single long time evolution when 
subsequent configurations are used in constructing the statistics. Our "statistical system" is now a single degree of 
freedom, the order parameter of the lattice system, interacting with all other (k ^ 0) modes. 

The effective equation of motion can be thought to result from the application of a "molecular dynamical rcnormal- 
ization group" , to our microscopical equations. The blocking in space is performed by projecting the field configuration 



^(x, i) on the OP $(t). It represents the infrared (IR) end point of such a blocking whose effective theory is now 
reconstructed from the actual time dependence found numerically. Combining ergodicity of the full system with the 
renormalization group (RG) concept wc arrive to the conclusion that ensemble averages of any OP-function coincide 
with time averages of the same function. 

In view of the RG concept we look for an effective equation of motion where the value of OP is determined exclusively 
by its values preceding in time. Assuming the absence of long memory effects a "gradient" expansion in time can be 
envisaged, leading to a local differential equation. We introduce into this phenomenological "Newton-type" equation 
a term violating time-reversal invariance. We are not able to derive it from the original system, we only wish to test 
its presence. The sustained motion of OP, however, requires, in this case the presence of a random "force" term, too. 
The deterministic part of the "force" is expected to be related to the equilibrium effective potential, since this object 
determines the stationary probability distribution for the OP near equilibrium. 

The above considerations lead us to write down the equation of motion, which represents a linear relation between 
the acceleration and the velocity of the order parameter: 

®d + Vd($d)$d ~h d + dV fj^ d ' = &, (4-1) 

where Cd is a noise term. In the corresponding dimcnsionless equation of motion the following new rescalcd quantities 
will appear: 




r) d = \m\r,, C d = C\ m \-\-. (4.2) 



The fitting procedure for the coefficients on the left hand side of Eq. ( |4.l| ) was the following. For a given interval 
of time It the region of the order parameter space visited by the system {<I>(i)|t £ I t } was divided into small bins. 
Having defined the time set T$ b = {f £ It\$b < $(t) < $6 + A$} corresponding to a given bin, the linear relation 
$(£{,) = —T)(^b)^(tb) — fi^b) was fitted with the method of least squares using the 3?(i&), $(£&) data measured at time 
moments if, belonging to T<$, b . We have obtained in this way the coefficient functions ^(^f,), /($&)• Once the functions 
r}($) and /($) = -h+V^($) are determined, we evaluate for each time i the expression $(t) + r]($(t))$(t) + /($(t)). 
Its actual value determines the random noise function £(i), whose statistical features (autocorrelation) should be 
extracted from the data. 

The effective force /($) calculated from the time-average of the oscillatory motion around the equilibrium, is ex- 
pected to agree with the force coming from the theoretically determined finite temperature effective potential calculated 
perturbatively in the cut-off two-dimensional field theory for some appropriately chosen value of the temperature |f6| . 
With one- loop accuracy the expected equality reads: 



f (f^*) measured h -\- 
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+ 0{T 2 ). (4.3) 

theory 



The expressions on the right hand side of this equation are connected to the dimcnsionfull quantities of the original 
one-loop computation in the following way: 

T = 7T^ T ^ v ^) = -\^ + W=aT^ v ^ A = -. (4.4) 
6|mp 2 4 6\m\ z a 

The dimcnsionless temperature is defined by the time-average of the kinetic energy based on the assumption that in 
the effective theory of the OP it has the usual expression in terms of $(i): 

km \ /Vi(*(f,x)) a = ip* = 5. (4.5) 



t^oo t J Q 2 v v " 2 2 

One should note that this definition in the dimcnsionfull version implies a "Boltzmann-constant" , |m| 2 multiplying 
T d . 

By measuring the order parameter average for different values of the external field h we obtain the magnetisation 
curve of the system. In the numerical work we restarted the computation at different values of the external magnetic 
field but one might as well change h adiabatically and measure the force law at each (quasi)cquilibrium point. The 
results should agree in the stable regime, up to the phenomenon of the hysteresis. The resulting curve can be viewed 
as the numerical Legendre transformation by identifying the external source h with the derivative of the effective 
potential, 

-h + V^ s (($) Kmeasured ) = 0. (4.6) 

It was found that /($) measured is always vanishing at $ = (Q)h, measured, thus the force acting on the order parameter 
at the equilibrium position induced by the external source h is indeed always — ((<!>);,)• 
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FIG. 3. The force /(<&) as the function of the order parameter on 40 x 40 lattices with T = 1/30 and 1/3. The force is 
measured by shifting the center of motion to different $ values by applying appropriate h fields to the system. The finite size 
effects are negligible, see text. For comparison we display the force arising from the derivative of the one-loop and for T — 1/3 
also of the two-loop expressions of the effective potential. 

The two sides of the relation ( |4.3|) are shown in Fig]|. To the values of the external field used in preparing Fig.|| 
(W6 = -0.5,-0.02,0,0.5,1,2) single runs were selected by the requirement that the system stayed in the metastablc 
vacuum up to the time 10 6 . The force shown in Fig.|| demonstrates that not only the equilibrium positions but also the 
fluctuations of the OP are governed by the static effective potential. Similar measurements, performed on 100 x 100 
lattices showed no finite size effects in the stable regime, h > 0. The error of /($) is not shown in the Figure, since 
the typical values (~ 0.004 for T = 1/3 and ~ 0.002 for T = 1/30) are too small to be displayed. 

Another piece in the effective equation of motion (4.1) is the friction term whose presence indicates the dynamical 
breakdown of the time inversion symmetry. The friction coefficient ri(<fr) proves clearly non-vanishing and shows 
only weak dependence on the actual value of the OP around its equilibrium position. The breakdown of the time- 
reflection symmetry in a closed system must arise only in presence of infinitely many degrees of freedom. Till then 
only statements on the Poincare-time can be made. Thus our non-zero results for r\ require further clarifications. 

The point is that there are two types of infinities, controlled by the temporal UV and the IR cutoffs, respectively. 
The spontaneous symmetry breaking is driven by the IR modes, and the non-trivial minima of the potential energy 
arise from the presence of infinitely many degrees of freedom in the IR (thermodynamical limit). On the contrary, the 
dynamical symmetry breaking |T7[ is the result of the effects of the derivative terms in the action and the infinitely 
many UV modes (continuum limit). The breakdown of the time inversion, being related to a time-derivative term in 
the effective equation of motion, should come from the UV, the short time behaviour of the system. 

In fact, one expects no friction when the UV cutoff, at is so small, that not enough energy can be dissipated during 
such a short time. In quantitative terms one should have 



a t > t 



2n 



0{T), 



(4.7) 



where r is the time scale of the fastest mode in the system. The right hand side relation of Eq.(4.7) gives an estimate 
of the maximal frequency from the free dispersion relation p\ = 4(sin 2 p x a/2 + sin 2 p y a/2)/a 2 + M 2 (T) of the lattice 
hamiltonian system for fixed a in the limit at — > 0. The fast modes absorb energy from the OP in a single time step 
for at > t, and the friction term should appear in the effective equation of motion. 
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FIG. 4. The friction r]({&)h=o) as the function of the timelike lattice spacing bt = na t . The time scale of the fastest mode 
is r w 2 and 0.22 for a = 1 and 0.1, respectively, according to the perturbation expansion. Typical error bars appear on the 
curves both for the IR and UV regimes. 



With help of the temporal blocking 
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one can construct the trajectories corresponding to larger values of the time cut-off, bt = nat. Such blocking 



was performed in time up to n = 2000 and a discrete time equation of motion of the form Eq.(4T) was reconstructed 
for the blocked trajectory <5 n (t). The non-trivial dependence of r\ on the temporal cutoff at is shown in Fig. ||. 

One can distinguish two regimes separated by a crossover apparently independent of the lattice size, located at 
b t ~ 0.2t. A scaling behavior is observed on the UV side, where ij(bt) tends to zero with a critical exponent close to 1 
and the force being constant. This result should be independent of the actual blocking details. In the other regime, on 
the IR side r\ goes over first into a ^-independent regime, corresponding to the saturation of the energy transfer from 
the OP and giving a stable, microscopic definition of the friction coefficient. Finally, in the far IR part a qualitatively 
different oscillatory behaviour sets in. 

The location of the crossover from the UV scaling regime to the plateau can be understood by writing the fluctuation 
dissipation theorem in the corresponding discretised form: (£ 2 ) = 2rjT/bt- The linearly increasing regime of r](bt) 
implies constant second moment for the noise. When 7/(6* ) reaches the plateau the second moment of the noise 
decreases like l/b t . The crossover therefore is located at the autocorrelation time scale of the noise. 

A qualitative interpretation of the oscillatory IR regime can be based on the observed small amplitude beating 
phenomenon in the OP trajectory. This can be recognized by closer naked eye inspection of the left side of Fig.[j], 
which persists further also on the right side. It is reflected in the OP autoc orre lation function, too, since in the course 
of the blocking an interference effect occurs on the right hand side of Eq^l^) due to this regularity. This feature is 



relevant to the value of 77(64), responsible for the decay of all kinds of fluctuations. 

The appareance of peaks in 77(64) at both the maximal destructive and constructive interferences can be modeled 
semi-quantitatively by identifying the beating part of the OP-motion Rci5$ = Rci5$o exp(iu>t) with the stationary 
solution of a single weakly damped driven harmonic oscillator, <5$ + r]8$ + cjq(5$ = /o exploit). The blocking acts on 
the trajectory £$(7;) as 5&(t) — > u6$(t), where u = (exp(iujbt) — l)/iu>bt. It does not change the relative phase of the 
driving force and of the blocked oscillation amplitude, leading to a relation between the parameters of the original and 
the blocked equation of motion: 
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(4.9) 



The friction coefficient for the blocked trajectory turns out to be 77 = (77 + culm(u 2 )) /TLe(u) . It is easy to see that 
Re(w) is vanishing at maximal constructive and destructive interferences. This provides singularities in 77. Whenever 
Re(w) = we have Im(w 2 ) = and numerator changes sign in the vicinity of the singularity. Thus 77 > apart for 
a short time interval around the singularities where the non-harmonic features should stabilise the fluctuations and 
keep 77 > 0, as observed in our simulation. 



The a dependence appearing in Fig.|] arises from the following two effects. One is that for larger a the maximal 
oscillation frequency is smaller and the time resolution of the system becomes cruder. Another is that larger a represent 
bigger physical volume, many more soft modes and less harmonic system, which tends to invalidate the simple picture 
based on a single harmonic mode. As a result the effects of oscillatory nature will be smeared, as one clearly recognizes 
in the figure. 



The quality of any proposed deterministic equation of motion (e.g. equations similar to Eq.(4.1) with zero on the 
right hand side), can be judged by the amplitude and the autocorrelation of its error term, the noise term £. The 
amplitu de o f the noise was found at least two-three order of m agni tude below the average level of the force as fitted 
to Eq. (4.1). The autocorrelation function of the noise of Eq. (4.1) appeared to be local, approximately of the form 
(Cd(t)Cd{t')) ~ S"(t — t'). The status of the fluctuation-dissipation theorem will be investigated for more complex field 
theoretical systems in future investigations. 



The sclfconsistcncy of the definition of the temperature in Eq.(4.5) and also the establishment of the thermal 
equilibrium can be tested further by plotting histograms for the following quantities: 
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Et = Ek + E p , 



(4.10) 



where /i 2 is the slope of the force as the function of the order parameter, determined numerically. Typical results 
for the energy-histograms are shown in Figj^. It shows perfect agreement of all slopes and good agreement with the 
expectations based on equipartition of the energy between the kinetic and the potential parts. We have checked for 
several temperatures, that the temperature determined by these histograms and from the |k|-spectra of the kinetic 
energy agree. This piece of information is parallel to the recent detailed investigations of the thermalisation in (1 + 1)- 
dimcnsional $ 4 -theories |Ts|J20|]. 
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FIG. 5. The histogram of the energies Et, Ek and E p of Eq.((4.1C)) on a 256 x 256 lattice with h = 2/\/6 



It is worthwhile noting that the relaxation time from a given initial condition to the thermally distributed state 
like in this figure takes at least one order of magnitude longer time for systems in the symmetric phase compared to 
the spontaneously broken case. At higher energy density (temperature) one would expect larger collision frequency, 
therefore shorter thermalisation time. The opposite result hints to the importance of slow, soft modes, whose presence 
is due to the symmetry breaking mechanism. We shall argue in the next Section that these modes are responsible for 
the realization of the Maxwell-cut in the potential term of the effective equation of motion for the OP. What we find 
remarkable is that these modes are present not only in the mixed phase but also near (meta)stable equilibria, among 
the dynamical fluctuations around the ordered vacuum in contrast to the massive perturbative excitation spectrum in 
the equilibrium. 

We complete the analysis of this section with the remark, that after the transition to the stable vacuum shortly a 
thermaliscd distribution is recovered for the OP at somewhat increased temperature which agrees with the value given 



by the equipartition. 



V. JUMP FROM THE METASTABLE TO THE STABLE VACUUM 



As it has been emphasized in the Introduction, there are (at least) two different descriptions of the transition 
from the metastable state to the stable one. In the first approach an expansion around the lower pass, the critical 
bubble (bounce) configuration yields a detailed space-time picture and the transition rate by means of the analytic 
continuation of the potential experienced by the OP near the (meta)stable position [|||. 

Another possibility is to provide for the OP a probabilistic description, obtained by the elimination of all other 
degrees of freedom in some kind of blocking procedure. This description is based usually on some underlying Master 
equation for the probability distribution of the OP and the resulting Fokker-Planck equation The transition 

to the stable state appears in this approach formally as a tunneling solution of the Fokker-Planck equation. The 
probabilistic feature of the dynamics of the OP is supposed to arise from assuming a statistical ensemble of the initial 
conditions. 

In our simulation wc find results analogous with the predictions of the probabilistic description by analysing the 
OP-motion starting from a single, well defined initial condition. One might wonder at this point if it is possible to 
understand the probabilistic tunneling of the OP by following the system from a unique initial condition. The self 
averaging in time can not be used for this argument since such a transition occurs only once during the evolution in 
a (quasi)irreversible manner. We have to develop a third approach to the metastable— >stable transition. 

In general, the effective equation of motion for (f>(t) reflects the typical landscape of the microscopic potential 
energy functional around the actual point $(x, t) in the configuration space. The classical origin of what appears as 
a tunneling on the Fokker-Planck level must be the arrival of $(x, t) to the vicinity of some narrow valley opening 
up towards the stable vacuum. In traversing this valley the landscape changes and the typical fluctuations will be 
different from those felt in the metastable regime. The constants parametrising the effective equation of motion must 
reflect this change. 

Our goal in this Section is to construct an effective description of the transition to the stable state by carefully 
tracing the time evolution of the OP. This will be achieved by projecting the microscopic equation of motion onto the 
homogeneous mode and phcnomenologically parametrising it similarly to Eq.(4.1). 

As long as the system is far from the narrow valley of the instability the force is time independent and agrees with 
the force derived from the perturbative effective potential according to the part of Fig. || corresponding to $ > 0.9. 
When the system arrives close to the entrance of the unstable valley, ($ « 0.8), the soft modes start to be important. 
This is reflected in the slight glitch in the leftmost piece of the measured force law in Fig. |3|. A sequence of glitches 
results in a situation depicted in Fig. ||. It shows the force acting on the OP in three successive time intervals preceding 
the event of tunneling for h = — 0.04/-V/6. The /($) curves determined using Eq. (4.1) in the disjoint time subintervals 
coincide within error bars for all, but the last one. In this last interval preceding directly the transition towards the 
direction of negative <& values, the fitted force bends down and its average becomes (a small positive) constant. This 
is characteristic feature of the instant when the system finds the entrance into the unstable potential valley. The 
vanishing of the force is an indication for the dynamical rea lization of the Maxwell-cut. The OP moves fast through 
the valley and the method of fitting the trajectory to Eq. ( |4. l[ ) for finding the force fails due to the insufficient statistics. 

The fluctuation moments depicted in Fig. |l| tell a bit more about this region. The increased values of the moments, 
the renormalized coupling constants in Wilsonian sense at vanishing momentum, indicate the enhanced importance 
of the soft interactions as the OP tunnels through the mean field potential barrier. This softening makes the OP 
fluctuating with larger amplitude. The valley of instability is in a surprising manner flatter than the typical landscape 
around equilibrium. The flatness along the motion of the OP (the mode k = in momentum space) comes from the 
Maxwell-cut. The average curvature of the potential in the transverse directions, i.e. for modes with |k| ^ 0, can 
be estimated by the second functional derivative of the two dimensional effective action with cutoff |k|, which can be 
taken as k 2 + V^'($). The increased MS fluctuation of the OP corresponds to a decrease in V^' =0 ($) in the valley. It 
pushes down V^'(<I>) also for small non-vanishing |k|, in the low momentum regime which is expected to be the most 
influenced by the changing landscape. 
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FIG. 6. The force /($) as the function of the order parameter on 100 x 100 lattice with h = — 0.04/v6- Each curve 
corresponds to time intervals of length 500|m| _1 . The interval corresponding to the continous curve ends 1000|m| _1 before 
the transition, the long dashed curve refers to the next time interval and the short dashed curve is the force measured during 
the last time interval preceding the transition. One clearly observes the bending down left from the central OP-value in the 
metastable regime. Typical error bars are shown on the leftmost points of the three curves. 

The lesson of Fig. ||is that the potential itself should be considered as a fluctuating quantity. This helps to translate 
into mo re q uantitative terms the above qualitative points. We discuss the projection of the microscopic equation of 
motion (|2.5|) onto the zero momentum sector, 



= $- $ + $ 3 4 

= * + 



V V 

3(p 2 $ + tp 3 - h 



(5.1) 



where the symbol ip n means the space average of (p n , (tp v = 0,$(i,x) = $(t) + ip(x,t)). The instant potential 
introduced in the second line contains a deterministic piece, which is the sum of the tree-level potential and of the 
slowly varying part the second and third moments (eg. (f nV )det)- This last feature clearly appears graphycally in 
Fig.|l] and a simple model based on a two-phase model of the transition period will be constructed below to account for 
it. The remaining oscillating pieces of the moments provide the probabilistic fluctuating contribution to the instant 
potential: 



The additive (Co) and the multiplicative (£i) noises are given by the differences 



(5.2) 



( (t) = ^ 3 (t) - (<^ ) deU 



(5.3) 



Note that no friction terms can be introduced in a natural way into the system (5.1-5.5) 



Therefore we face the 
— r* 

r) appears 



intriguing question, how irreversibility is realised in such a system. The time-correlation matrix Q(t)Q (t 
to us to be the key object for its investigation, to which we plan to return in the future. 

In order to gain more insight how this works we build into Eq.(5.1) the consequences of the mixed two-phase picture 
of the phase transformation. The microscopical basis for this picture is provided by the thermal nucleation whose 
quantitative discussion is given for completeness in the Appendix. 

More specifically, we assume that the space can be splittcd into sharp domains (neglecting the thickness of the walls 
in between), where the field is the sum of the constant background values &o± and the fluctuations ip± around it, 



$±(x,i) = $ 0± + <p±(x,i). 



(5.4) 



We assume local equilibrium in both phases, based on the smooth evolution of the temperature as displayed in Fig. 



The actual value of the order parameter is determined by the surface ratio p(t) occupied by the stable phase: 



$(t)=p(t)$o- + (l-K*))*o+. (5-5) 



Simple calculation then yields 
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fo^WC*) = ° + - U ($ 3 _(x,*) - (x,t) V ) + *%{x,t) V - 3*(t)(^ V ) (fat (t) - <& 3 (t), (5.6) 

where the volume averages should be read off the corresponding equilibria on the two sides of the transition. If one takes 

the values of <J>o± , ip± from the respective equilibria determined in the same simulation, a quite accurate description 

of the shape of the two fluctuation moments arises in the whole transition region and its close neighbourhood using 

v 

<£>(£) to parametrize their t-dependence. (Note that $±(x,t) does not depend on time.) The deterministic part of 
the moments are therefore well-defined functions of $. In this way, a simple explicit construction can be given for the 
effective noisy equation of the OP-motion, if the above pieces arc supplemented by the correlation characteristics of 
the two kinds of noises. 

The RMS (root mean square) fluctuations found numerically on the transition part of the trajectory and its close 
neighbourhood are the following: 



\jQ{t) = 0.0063(10), = 0.0054(10), for N = 64 (5.7) 

(its magnitude increases with the size of the system) . The magnitude of their equilibrium cross correlation was found 
- 10" 5 . 

As a corollary of this construction one can demonstrate the absence of the deterministic part of the acceleration of 
the order paramete r in the trans ition period. 



Substituting Eq.(5.6) into Eq.(5.1), one finds for the deterministic part of the force, 

, s / $i(x,t) V - $ 3 (x, t) V \ $ 0+ $ 3 (x,£) y - $ _$i(x, t) 1 , . 

/w =(- i+ — )* + (5 - 8 > 

The average of the equations of motion in the respective equilibria, 

(*!(x,t)>-$o±-h = (5.9) 
implies the vanishing of the deter mini stic force in Eq. (|5 .8|) , when exploiting the equality of the volume and the 



ensemble average in this case. Eq. (5.1) transforms into 

*(t) + Co(*) + 3&(t)$(t)=0. (5.10) 

This is the dynamical realization of the Maxwell-construction holding when the mixed phase model with local equi- 
librium is valid. 



If the force is calculated from the full instant potential, ^( c &) ms t m Eq- (5-2), it depends parametrically on the 
— v — V 

moments <p 2 (t) and (p 3 (t). The approximate trajectory ( l ) i ns t(^); defined by minimising V"($) ms .t. with respect to 
where the moments are taken from the numerically determined time evolution, reproduced accurately the observed 
OP, with the following RMS/unit time: 

VSrS? =0.0014(5), for7V = 64. (5.11) 

This construction intcrpretes the OP-trajectory as a continous deformation of the instant potential with the OP 
" sitting" permanently in the actual minimum. Notice that such motion is possible only if the OP continously undergoes 
some sort of dissipation. 

The vanishing of the acceleration was checked by comparing the computer generated OP trajectory in the transition 
region with a ballistic motion in a viscous medium. In particular it was tested that the ratio h/<& is nearly time and 
h independent for small enough h. The dynamical friction measured by the above ratio tends to a constant with 
decreasing h at fixed lattice size. Although the error for its value in each individual transition is rather small, the 
central values obtained in different runs fluctuate quite strongly, which leads eventually to a large error of the mean 
calculated as an average of runs with different initial configurations having the same energy density. 

It is worth to note that the order of magnitude of these "renormalized" , i.e. IR determined values of the friction 
coefficient agree with the peak value in Fig. |i| for N = 64. The "running" friction coefficient determined by the 
blocking in time, however, drops as bt is increased and takes extremely small values at bt ~ 20, the average time length 
of the ballistic fits for the transition. This serves an example for the dependence of the renormalization group flow 
on the details of the blocking. It remains to be understood whether the agreement between the peak value and the 
ballistic fit is an accident or follows from the internal dynamics. 



VI. CONCLUSIONS AND FUTURE DIRECTIONS 



In this paper we have investigated in detail the decay of the false vacuum in a classical lattice field theory based 
exclusively on the effective theory of the order parameter. Two versions of the effective theory were reconstructed 
from its trajectory derived from the microscopical equations of the theory. The first refers to the (meta)stable branch 
of the motion. The second one which takes into account the existence of a mixed phase during the transition period 
describes very well the transition together with its neighbourhood. The first equation has the form of conventional 
mechanical motion taking place in a dissipative noisy environment. The dissipation, the dynamical breakdown of the 
time inversion symmetry was found only for times longer then the minimal microscopic time scale of the system, the 
autocorrelation time of the noise. 

As a corollary we find that the effective phenomenological OP-theory for the decay of the false vacuum through a 
narrow but flat valley demonstrates the presence of a dynamical Maxwell-construction. 

The next possible directions of extending this work include the investigation of the effect of the quantum fluctuations 
on the large amplitude thermalization of the one-component $ 4 theory ]2l| , p2| ] . The effective equation of motion for 
the order parameter of a quantum field theory can be studied by breaking up the field $(x, t) into a "classical" and a 
quantum part: 

$(x,t) = * cl (x,t) + 0(x,t), (*(x,i)> = $ d (x,t). (6.1) 



The quantum effects enter the equation of $ c ](x, t) through (</> n (x, t)), n = 2, 3 analogously to Eq. fl5.l| ). The method 
of the mode-function expansion |2l|-|23]| displays explicitly the way the quantum effects (~ Ti) enter the "classical" 
equations of motion. In this method one solves the time evolution of the coefficient functions of the expansion 

0(x, t)=J2 [/q(x, t)oQ + /£(x, t)a+\ , (6.2) 
Q 

where [<xq, Oq,] = <5qq' . The initial conditions are fixed by requiring the canonical commutation relations to be fulfilled 
at t = 0: 



The quantum effects relative to Eq. (2.5) show up as 0{%) corrections to the coefficients of the effective equation of 

v 

motion of "^(x, t). The adequate form of the effective equation of motion for '^(x, t) represents a further challenge 
ll- 

The extension of our results to classical field theories with continous internal symmetry describing the dynamics of 
the inflaton field coupled to the Higgs+Gauge system P5p9] is also of actual interest. 



APPENDIX: THE NUCLEATION PICTURE 



In this section we analyse the transition using the more conventional statistical approach of thermal nucleation 
theory. 

From the study (on lattices up to N = 512) of the detailed microscopical field configurations it turned out that the 
phase transformation starts by the nucleation of a single bubble of the stable ($ < 0) phase. Late-coming further large 
bubbles are aggregated to it as well as the small size (consisting of Ububble < 50 — 60 joint sites) bubbles. Following 
the nucleation the expansion rate of the large bubble governs the rate of change of the order parameter. To a very 
good approximation each individual transition could have been characterised by a constant value of $ in this interval. 
Two mechanisms are known to lead to this behavior. The first is scattering of hard waves ("particles") off the bubble 
wall, while in the second the expansion velocity is limited by the diffusive aggregation of smaller bubbles |27[ ]. The 
latter process seem to be the dominant in the kinetic Ising model fl5| . 

The statistics of the "release" time of the supercritical bubble from the metastable state shows at first sight rather 
peculiar characteristics. The binned histogram for larger values of h shows an asymmetric peaked structure which 
apparently deviates from the exponential distribution characterising the thermal nucleation scenario (Figj^) . The very 
early transitions (t < t rnax ) seem to be suppressed for small values of h. We expect they correspond to transitions, 
which happen before the system reaches the metastable state. 

The fall-off for t > t max starts nearly exponentially, but the histogram develops a very long tail when h — ► 0. The 
bigger sample is used for estimating the probability distribution, the more suppressed is the weight of these events 
in the normalised distribution. In practice this means a longer time interval where a good linear fit can be obtained 
to the log-linear histogram. Eventually, the separation of a clean exponential signal is possible, and the slope can be 
compared with predictions of the nucleation theory. 
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FIG. 7. The normalised "release" time histogram counting the frequency of a certain time moment-bin when the order 
parameter first takes a negative value. The figure represents the statistics based on 750 events at h = —O.OS/VG on an TV — 64 
lattice. 

The nucleation rate is proportional to the volume of the system (^V 2 ). In Fig.|| we see evidence for the size 
independence of the transition rate per unit surface. The logarithm of this quantity can be estimated following the 
standard nucleation theory |^,^,^8| : 



lnT = K 



T 



(6.4) 



S-2 is determined approximately by the action of the bounce solution, connecting the two vacua. In this investigation 
we have studied at fixed temperature the /i-dependence of S%/T, assuming the ^.-independence of K. In the thin wall 
approximation the following expression is used based on the tree level potential: 
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(6.5) 



(the first term on the right hand side corresponds to the volume energy of a bubble of radius f, the second one to its 
surface energy). For the critical bubble size a very simple expression is obtained for the exponent of the rate in terms 
of dimcnsionlcss quantities: 



5 2 (thin wall) 16tt \m\ 5 1 
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(6.6) 



The linear dependence on l/h is fulfilled in our numerical calculations very well, but the predicted action is more 
than one order of magnitude larger than what can be derived from the slope of Fig.||: S2/T (measured) « 0.1/(/iT). 
If one relaxes the thin wall approximation and solves numerically the two-dimensional bounce equation directly for 
several h, one finds S 2 (bounce)/T = 0.78/ (hT). 
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FIG. 8. Nucleation rate for unit lattice surface versus the reciprocal of the magnetic field for different lattice sizes, 
error bars are smaller than the size of the symbols representing the data points. 



The 



Further improvement can be obtained by applying the temperature corrected effective potential in the bounce 
equation. We have determined the parameters of the T-dependent potential directly from our numerical calculation in 
the following way. The restoring force was measured for a certain h around both the stable and the mctastablc minima 



as decribed in section 4. Next an interpolating fit has been constructed of the form: h e Q - 



TO eff $ + A 3$ 2 



It turned out that in the best fit A3 « is fulfilled always, while the values of the other coefficients are not too 
far from their tree level values. Then a bounce solution can be built on the corresponding (real!) potential which 
includes the temperature corrections. The final result is quite close to the measured value of the rate logarithm: 
S2 (T, bounce)/T = 0.29/(/iT). The same real interpolation can be built on the neighbourhood of the minima of the 
2-loop T-dcpcndcnt effective potential, leading to ^(T, bouncc)/T = 0.2/ (hT). By common experience in the surface 
tension simulations a factor of 2-3 difference in £2 is expected to arise relative to the mean field theory. 

We conclude, that the finite temperature corrections are important for quantitative treatment of the nucleation rate 
within the thermal nucleation theory. 

Next we turn to the discussion of the possible nucleation threshold at small h. The average "release time" increases 
for fixed lattice size with decreasing h (Fig.^J). On a lattice of fixed size we found small h values for which we could 
not detect any transition, what makes very probable the existence of a threshold value of the external field h t h(N). 
We did not attempt to locate this value beyond the simple hyperbolic fits to the few largest (t re i ease ) values. The 
value of hth(N) remains stable when the smallest h is excluded from the fit, therefore we conclude that a threshold 
magnetic field exists similar to the case of the kinetic Ising model ]l5| |. The value of h t h decreases when the lattice 
size is increased. Intuitively we expect h(oo) = 0, but with the three lattices studied by us this conjecture cannot be 
demonstrated. 
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FIG. 9. The average "release" time as a function of h measured on N = 64, 128, 256 lattices. 



ACKNOWLEDGEMENTS 



The authors are grateful for valuable discussions to C. Gagnc, J. Hajdu, A. Jakovac, Z. Racz and H.J. de Vega. 
They acknowledge the use of computing resources generously provided by the University of Bielefeld and the Dept. of 
Computer Science and the Inst, for Theoretical Physics of the Eotvos University. This research has been supported 
by CNRS and the Hungarian Science Fund (OTKA). 



[1] Nucleation, eds. A. Z. Zettlrmoyer, M. Dekker (Academic Press. N.Y. 1969); F. F. Abraham, Homogeneous Nucleation 
Theory, (Academic Press, N. Y. 1974); J. D. Gunton, M. San Miguel, Paramdeep Sahni, in Phase Transitions and Critical 
Phenomena, vol. 8. eds. C. Domb., J. L. Lebowitz (Academic P ress, N. Y., 1983) 

[2] for a recent introduction, see D. Boyanovsky and H. de Vega, |iep-ph/9909373 , to appear in the Procs. of the IV Paris 
Colloquium on Cosmology 

[3] K. Binder, D. Stauffer, Adv. Phys. 25 (1976) 343 

[4] J. S. Langer, Ann. Phys. (NY) 41 (1967) 108 

[5] S. Coleman, Phys. Rev. D15 (1977) 2929, S. Coleman and F. Do Lucia, Phys. Rev. D21 (1980) 3305 
[6] J. Alexandre, V. Branchina, J. Polonyi, Phys. Lett. B445 (1999) 351 

[7] D. J. E. Callaway and D. J. Maloof Phys. Rev. D27 (1983) 406, L. ORaifaertaigh, A. Wipf and H. Yoneyama, Nucl. Phys. 
B271 (1986) 653, M. Hindmarsh, D. Johnston, J. Math. Physics A19 (1986) 141, V. Branchina, P. Castorina, D. Zappala, 
Phys. Rev. D41 (1990) 1948, A. Ringwald, C Wetterich, Nucl. Phys. B334 (1990) 506 K. Cahil Phys. Rev. D52 (1995) 
4704 

[8] D. Boyanovsky, H.J. de Vega, R. Holman and J. Salgado, Phys. Rev. D59 (1999) 125009 

[9] M. Alford and M. Gleiser, Phys. Rev. D48 (1993) 2838 
[10] M. Gleiser and R.O. Ramos, Phys. Rev. D50 (1994) 2441 
[11] M. Gleiser, Phys. Rev. Lett. 73 (1994) 3495 

[12] A. Berera, M. Gleiser and R.O. Ramos, Phys. Rev. D58 (1998) 123508 
[13] E. Brezin and J. Zinn- Justin, Nucl. Phys. B257 (1985) 867 
[14] K. Binder, Z. f. Physik B43 (1981) 119 

[15] K. Binder and H. Muller-Krumbhaar, Phys. Rev. B9 (1974) 2328 

[16] J.I. Kapusta, Finite Temperature Field Theory, (Cambridge University Press, Cambridge, 1989) 

[17] V. A. Miransky, Phys. Lett. 91B (1980) 421, J. Polonyi, in the Proceedings of the International School of Physics, Enrico 

Fermi, Course CXXX, on Selected Topics in Nonperturbative QCD, Soc. It. di Fisica 
[18] N. G. van Kampen, Stochastic Processes in Physics and Chemistry, North Holland, 1981, see also A. Linde, Nucl Phys. 

B372 (1991) 421 



[19] G. Parisi, Europhys. Lett. 40 (1997) 357 

[20] G. Aarts, G.F. Bonini and C. Wetterich HD-THEP-00-20, frep-ph/0003262 

[21] D. Boyanovsky, H.J. de Vega, R. Holman, D.-S. Lee and A. Singh, Phys. Rev. D51 (1995) 4419 

[22] S. Khlebnikov and I.I. Tkachev, Phys. Rev. Lett. 77 (1996) 219, S. Khlebnikov, in Procs. of the Int. Conf. "Strong and 

Electroweak Matter '97", eds. F. Csikor and Z. Fodor, World Scientific, Singapore, 1998, pp. 69-88 
[23] J. Baacke, K. Heitman and C. Paetzold, Phys. Rev. D57 (1998) 6398 
[24] Sz. Borsanyi, A. Patkos, Zs. Szep work in progress 

[25] J. Garcia-Bellido, D. Grigoriev, M. Kusenko and M. Shaposhnikov, Phy s. Rev. D60 (19 99) 123504 

[26] D. Boyanovsky, D. Cormier, H.J. de Vega, R. Holman and S.P. Kumar, hep-ph/9801453] 

[27] T.A. Witten and L.M. Sander, Phys. Rev. Lett. 47 (1981) 1400 

[28] M. Gleiser, G. C. Marques and R. O. Ramos, Phys. Rev. D48 (1993) 1571 



